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Abstract 

We re-evaluate the nuclear suppression of forward 7/t/t production at high energy in the Color Glass Condensate 
framework. We use the collinear approximation for the projectile proton probed at large x and an up to date dipole 
cross section fitted to HERA data to describe the target in proton-proton collisions. We show that using the Glauber 
approach to generalize the proton dipole cross section to the case of a nucleus target leads to a nuclear modification 
factor much closer to LHC data than previous estimates using the same framework. 
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1. Introduction 

The study of forward Jjtf/ production in high energy 
proton-proton and proton-nucleus collisions, which 
probes the target at very small x, can provide valu¬ 
able information on gluon saturation. Indeed, the charm 
quark mass should be small enough to be sensitive to the 
saturation scale. On the other hand, it is large enough to 
provide a hard scale and thus to allow the use of a weak 
coupling treatment. It also has a clean experimental sig¬ 
nature so this process has been the subject of many ex¬ 
perimental studies to date. 

In this work we will study 7/i/f production at for¬ 
ward rapidities in proton-proton and proton-nucleus col¬ 
lisions at the LHC in the Color Glass Condensate (CGC) 
framework, using the color evaporation model (CEM) 
to treat hadronization. Since we work at forward rapid¬ 
ity, where the projectile is probed at large x and the tar¬ 
get at small x, we will use the “hybrid model” in which 
the projectile proton is treated as dilute and is described 
in terms of usual collinear parton distribution functions 
(PDEs). The collinear gluon emitted can then split into 
a cc pair either before or after the interaction with the 
target. These partons are then assumed to eikonally 
interact with the target, picking up a Wilson line fac¬ 


tor in either the adjoint or the fundamental representa¬ 
tion, depending on the particle. The cross section for 
cc pair production is then described in terms of Wil¬ 
son line correlators containing the information on the 
dense target. The same Wilson line correlators appear in 
calculations of other processes, such as total DIS cross 
sections, single and double inclusive particle production 
in proton-proton and proton-nucleus collisions, diffrac¬ 
tive DIS and the initial state for hydrodynamical model¬ 
ing of heavy ion collisions. This framework has thus a 
broad range of applications. 

The modification of 7/i/^ production cross section 
in proton-nucleus compared to proton-proton collisions 
has been previously studied in the CGC framework |[T|. 
However it was found that the nuclear suppression pre¬ 
dicted by this calculation was much stronger than mea¬ 
sured later at the LHC. In this work we re-evaluate this 
quantity in the same collinear “hybrid” framework, us¬ 
ing a more careful treatment of nuclear geometry neces¬ 
sary to go from the description of a proton to the one of 
a nucleus in the CGC framework. This is motivated by 
the fact that it was observed for example in single inclu¬ 
sive light hadron production Q that the disagreement 
of previous CGC calculations lO with LHC data was 
mostly due to nuclear geometry effects. We also use the 
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more recent dipole cross sections which were obtained 
in Ref. E). 


2. Formalism 


In this work we use the simple color evaporation 
model (CEM) to describe the hadronization of the cc 
pair into a 7/i/f meson. We note that it is also possible 
to treat hadronization in a more elaborate way, for ex¬ 
ample by using an expansion in terms of non-relativistic 
QCD as was done in Ref. ID, but here we focus on the 
importance of nuclear geometry. In the CEM a fixed 
fraction of the cc pairs, produced either in the color sin¬ 
glet or octet state, whose invariant mass is below the 
D-meson threshold, is assumed to hadronize into Jl\}/ 
mesons. The differential cross section with respect to 
the transverse momentum and the rapidity Y of the 
produced 7/i/r reads 


d^PxdT ]Am) d2p_LdFdM2’ 


( 1 ) 


Where nic is the charm quark mass, nio - 1.864 GeV 
is the D meson mass and ^ip'^ydM- cross section 
for cc pair production with transverse momentum 
rapidity Y and invariant mass M. The nonperturbative 
constant Fjj^ in Eq. [T]is related to the probability for 
a cc pair to transition to a 7/i/r. In this work we will be 
mostly interested in the nuclear modification factor RpA, 
defined as 


1 dcr/d2p^dy| 

= T(2) 

A do-/ d2PA dT 

ipp 

where do-/ d^Pj^ dTl and do-/ d^Pj^ dT| , are the 
cross sections in proton-proton and proton-nucleus col¬ 
lisions respectively. Therefore for this observable 
plays no role and we don’t need to fix it. 

The formalism for gluon and quark pair production 
in the dilute-dense limit of the CGC has been studied 
in detail in Refs. Elia (see also Ref. Q) and used in 
several works, such as EEiiiiiia. In this framework, 
when using the collinear approximation to describe the 
gluon emitted by the projectile, the cross section for cc 
pair production reads, in the large-A/; limit IfTI : 

do-cg _ alN^ 1 r Scoii(/7y -I- 

d^Px d2^^ Ayp (22r)2 j {p^ + 

X i _L(P± + x\Gp{xx , Q^) , (3) 

yt X2 


where and are the transverse momenta of the 
quarks, yp and their rapidities, = J 

and dp, = - 1 is the dimension of the adjoint repre¬ 

sentation of SU(A/;). The expression for the “hard ma¬ 
trix element” Hcoii can be found in Ref. m. The longitu¬ 
dinal momentum fractions probed in the projectile and 
the target, x\ and X 2 , are 



(4) 


The propagation of the cc pair in the color field of the 
target is described by the function 

= f S,(kJ S,(L-kJ , (5) 

where b± is the impact parameter. The function Sy(kj_) 
is the fundamental representation dipole correlator in 
the color field of the target and it contains all the in¬ 
formation about the target. It reads 

Sy(kJ^ J dVye'*--^-5,(ry) , (6) 

with 

Sy(x^-yJ^ ^(TiUHxJUiyj), (7) 

-/Vc 

where t/(jCj_) is a fundamental representation Wilson 
line in the color field of the target. 

In the case of a proton target, where there is no ex¬ 
plicit dependence of the dimensionless dipole amplitude 
on the impact parameter, the following replacement is 
made: 

f ^ Y , (8) 

where ctqII corresponds to the transverse area of the 
proton measured in DIS experiments. The function 
then reads in this case 

^ ^ Sy(kj Sy(h - kj. (9) 

^ L 4as 


Eor the description of the gluon distribution in the 
projectile proton Gp(xi,Q^), treated in the collinear 
approximation, we use the MSTW 2008 IfTTl LO 
parametrization since the remaining of our calculation 
is done at leading order. 


3. Dipole correlator 

Eor the initial condition to the running coupling 
Balitsky-Kovchegov equation lfT2l [T3l [T4l governing 
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the rapidity evolution of Sy(rj_) we use the MV® 
parametrization from Ref. m at initial rapidity Y - 
ln(l/xo) with xo — 0.01: 

5y=ini(^-±) = exp 

^0 

( 10 ) 

and the following expression for the running coupling: 

127T 

Q:s( 0 = --r- (11) 

(33-2iVy)log(^) 

The parameters = 0.060 GeV^, - 7.2, e^. = 18.9 
and cro/2 = 16.36 mb are obtained by fitting HERA 
DIS data 113 for < 50 GeV^ and x < 0.01. One ad¬ 
vantage of this parametrization over the AAMQS llT3 
one, similar to the one used in Ref. U, is that Sy(kj_) 
is positive definite at any rapidity, which is necessary to 
interpret it as an unintegrated gluon distribution. Never¬ 
theless we note that at LHC energies for which we will 
show results here, the dependence on the precise form 
of the initial condition is rather weak. 

While the dipole correlator for a proton target can be 
obtained by hts to DIS data, the extension to a nucleus 
needs to be based on some assumptions. In Ref. |[T] the 
dipole correlator for nucleus was obtained by using the 
same initial condition as in the proton case, but with an 
initial saturation scale scaled by a factor ~ Here 

we instead use the Glauber approach, in which the initial 
condition is given by 


rx^Q: 


sO 


In 


1 


kxIAgcD 


■+er 


SY^i^±(r±,bJ = exp 
-'0 


AT^bJ 


o-Q 

2 4 


X In 


1 

kxIAqcD 



( 12 ) 


where the transverse thickness function Ta is given by 
the standard Woods-Saxon distribution 


TAibJ 


r dz — 

1 + exp 


n 

d 


(13) 


with d = 0.54fm and Ra = (1.12A'/3 - 
Here n is a normalization factor fixed so that the integral 
of Ta over the impact parameter is unity. To compute 
observables with nuclear targets we use the expressions 
given in the previous section but instead of making the 
replacement we integrate explicitly over the impact 
parameter. An issue in this approach is that the nucleus 
becomes dilute close to its edge, with a saturation scale 
being smaller than the one of the proton. In this region. 


BpPb 



Figure I: Nuclear modification factor for Jlijj production as a function 
of Y. Data from Refs. I18II19I . 


which has a small contribution to the total cross section, 
we use the proton-proton result scaled so that - 1 ■ 
In this treatment the only additional input when going 
from a proton to nucleus target is thus the standard nu¬ 
clear density Ta- 


4. Results 

We now turn to our results. We will here focus on the 
nuclear modification factor R^a, for results on the cross 
sections for proton-proton and proton-nucleus collisions 
we refer to Ref. ini , where it was observed that the gen¬ 
eral shapes for the cross sections as a function of Ry and 
Y are in quite good agreement with experimental data, 
although they are affected by a quite large normalization 
uncertainty. These uncertainties cancel to a large extend 
in RpA, as can be observed from Figs.[^and|^where we 
show Rppb at yjsAiN - 5 TeV as a function of Y and Ry 
respectively. The uncertainty band in both figures corre¬ 
sponds to the variation of the charm quark mass between 
1.2 and 1.5 GeV and of the fact orization scale between 
Mj_/2 and 2Mj_ with Mj_ = + P\ and M is the 

cc pair’s invariant mass. The comparison with measure¬ 
ments at the LHC by the ALICE and LHCb collabo¬ 
rations shows that the values obtained here are much 
closer to the data than the previous work in Ref. fT], 
but they are still slightly too small to describe the data, 
especially at low Ry as can be seen from Eig. 


5. Conclusions 

In this work we have re-evaluated the CGC predic¬ 
tions for the production of 7/i/f mesons in high energy 
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^pPb 



Figure 2: Nuclear modification factor for 7/^ production as a function 
of Px (2 < F < 3.5). Data from Ref. (T^. 

proton-proton and proton-nucleus collisions, in partic¬ 
ular focusing on the treatment of the nuclear geome¬ 
try when going from a proton target to a nuclear one. 
The dipole cross section in the proton case is fully con¬ 
strained by DIS fits, while the only additional input 
used for a nucleus target in this approach is the standard 
Woods-Saxon distribution. Our results in this Glauber 
approach for the nuclear modihcation factor, which has 
a quite small uncertainty, are closer to experimental data 
than previous estimates. They nevertheless slightly un¬ 
derestimate this ratio, especially at low transverse mo¬ 
mentum. 
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